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Abstract. The relative orientation of clusters' major elongation axes and clusters' angular momenta is studied 
using a large N-body simulation in a box of 500 /i _1 Mpc base length for a standard ACDM model. Employing 
the technique of mark correlation functions, we successfully separated the correlations in the orientation from 
the well known clustering signal traced by the two-point correlation function. The correlations in the orientation 
are highly significant for our sample of 3000 clusters. We found an alignment of neighboring clusters, i.e. an 
enhanced probability of the major elongation axes of neighboring cluster pairs to be in parallel with each other. 
At 10 /i -1 Mpc separation the amplitude of this signal is ~ 10% above the value expected from random orientations, 
and it vanishes on scales larger than 15 ft~ 1 Mpc. The "filamentary" alignment between cluster's major elongation 
axes and the lines pointing towards neighboring clusters shows even stronger deviations from random orientation, 
which can be detected out to scales of 100 ft _1 Mpc, both in 2D and 3D analyses. Similarly, strong correlations of 
the angular momentum were seen. Also a clear signal in the scalar correlation of the absolute value of the angular 
momentum, the spin parameter and the mass was found. They extend up to 50 /i _1 Mpc and have an amplitude 
of 40%, 15%, and 10% above a random distribution at 10 /i~ 1 Mpc separation, respectively. 

Key words, large-scale structure of the Universe - methods: statistical - galaxies:clusters,general 



1. Introduction 

The study of orientation effects between galaxy clusters 
has a long and controversial history in cosmology. In a 
seminal study Binggeli ( 1982 ) claimed that galaxy clusters 
are highly eccentric and oriented relative to neighboring 
clusters if lying at separations smaller than 15 /i _1 Mpc. 
Further he found anisotropics in the cluster distribution 
on scales up to 50 /i _1 Mpc. Following studies found no 
or weak statistical significance for orientation effects be- 
tween neighboring clusters or between cluster orienta- 
tion and the orientation of the central dominant galaxy, 
Strublc fe Peeblesj fll985| ), fUn] ( |1987| ), and |Rhee fc 



cp 



Katgert (1987). Remarkable was the apparent absence 
(Ulmer et al. 1989) or weakness (jRhcc fc Latour 1991) 



new rich cluster catalogues as the optical Enacs survey 
(Katgert et al. 199£) and the X-ray based Reflex survey 
(Bohringer et al. 1998), the question of orientation effects 
in clusters should attract renewed attention. Sufficiently 
large and well defined cluster samples showing only weak 
contamination by projection effects seem to be necessary 
to clarify this uncertain situation. 

Strong stimulus to study orientation effects in clusters 
came from early ideas that a possible relative orientation 
between neighboring clusters or of clusters in the same 
supercluster should reflect the underlying structure for- 
mation mechanism. Binncy & Silk ( 1979| ) proposed that 
tidal interactions of evolving protocluster systems may 
lead to the growth of anisotropics of clusters and to rel- 



ative orientation effects. Later, van Haarlem ct al. (1997) 



of orientation effects in projected X-ray contours of clus- 
ters, but it should be noted that the cluster samples at 
this time were small. Analyzing a large set of 637 Abcll 



used numerical simulations of CDM models to demon- 
strate that clusters are elongated along the incoming di- 



rection of the last major merger. In the same spirit, West 



clusters, Plionie (1994) found highly significant alignment 
effects on scales below 10 ft, _1 Mpc that become weaker but 
extend up to 150 /i _1 Mpc. More objectively selected, but 
smaller cluster sample s seemed to put in to qu estion the re- 
ality o f this signa l, cp. Fong et al. ( 1990|) and Martin et al 
(|1995|) . However, |Chambers et~ai] (|2000| ) found significant 



nearest neighbor alignment of cluster X-ray isophotes us- 
ing data from Einstein and ROSAT. With the advent of 



(1994) found that clusters grow by accretion and merging 
of surrounding matter that falls into the deep cluster po- 
tential wells along sheet-like and filamentary high density 
regions. Therefore, the cluster formation is tightly con- 
nected with the supercluster network that characterizes 
the large-scale matter distribution in the universe. High- 
resolution simulations showi ng this effect are des c ribed by 
the Virgo collaboration, cp. Colberg et al. ( 2000 ). Onuora 
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& Thomas ( 2000| ) found a significant alignment signal up 
to scales of 30 ft, _1 Mpc for a ACDM model, whereas in 
a tCDM model the signal extended only up to scales of 
lSft^Mpc. 

To quantify the alignment of the galaxy clusters, we 
use a large ACDM simulation in a box of 500 ft, Mpc 
side length. We identify a set of 3000 clusters. As statisti- 
cal tools we employ mark correlation functions (MCF), as 
introduced to cosmology by Beisbart & Kerscher ( 2000) ). 
In this article we will extend this formalism to allow for 
vector valued marks. The direction of the major axis of 
the mass ellipsoid serves as the vector mark. Tightly con- 
nected with the elongation of clusters is its internal rota- 



tion. According to |Doroshkcvich| ( |1973|) and |Whitej (|198J), 
the primary angular momentum of bound objects is due 
to tidal interaction between the elongated protostructures 
after decoupling from cosmic expansion and before turn- 
around. More recent studies find that the angular momen- 
tum of dark matter halos is later modified by the merging 
hist ory of their bui l ding b locks, cp. Vitvitska et al. ( 2001 ) 
and Porciani et al. ( 2002a| ,|bl)- Therefore, we utilise the an- 
gular momentum as an additional mark for the study of 
the correlation of inner properties of simulated clusters, 



and wd compare it with the orientation effects. 



The plan ot the paper is as follows. In the next sec- 
tion, we describe our numerical simulation, the selection 
of a cluster sample and the precision with which we can 
determine structure parameters from it. Next we discuss 
the MCFs that are relevant for our studies. In particular, 
we use special MCFs for vector marks to quantify correla- 
tions of orientation. In Sect. |^ we investigate correlations 
in the spatial orientation of clusters both in 3D and in 
the projected mass distribution. In Sect. ^| we present a 
MCF analysis using the angular momentum, mass and 
spin taken as vector and scalar marks, respectively. We 
conclude with a summary of the results. 

2. Cluster sample in numerical simulations 

We utilize the AP3M code of |Couchmaij ( |l99lD to follow 
the dynamics of 256 3 particles in a box of 500 h~ 1 Mpc 
with periodic boundary conditions. We employ a cold dark 
matter model with a cosmological constant ft\ — 0.7, 
a matter density f2 m = 1 — CI a, and a Hubble constant 
H = 100ft, km/s/Mpc with ft = 0.7. The age of the 
universe in this model is « 13.5 Gyrs. The normaliza- 
tion, given by the linear mass variance of dark matter on 
8 ft, _1 Mpc scale, <j$ — 0.87, is in accordance with the four 
year COBE DMR observations as well as the observed 
abundance of galaxy clusters. The code uses glass-like ini- 
tial conditions, cp. Knebe & Mullei (1999). The initial 



power spectrum was calculated with the CMBFAST code 
( [Scljak fc Zaldarriaga 1996). We start the simulations at 



an initial redshift z = 25. Up to this time the Zeldovich 
approximation provides accurate results on the scales con- 
sidered here. We employ a comoving softening length of 
100 hr l kpc, and 1000 integration steps that are enough 
to avoid strong gravitational scattering effects on small 



scales, cp. Kncbc et al. (2000). With this softening length, 
the inner cores or break radii of the cluster sized halos are 
resolved. For our statistical analysis the big simulation 
volume provides us with a sufficient number of clusters 
embedded into the large scale structure network. In par- 
ticular an accurate representation of the large scale tidal 
field is important for getting stable orientation results. 
The particle mass in the simulation is 6.2 x 10 1 1 Ii~ 1 Mq, 
comparable to the mass of one single galaxy. For reliable 
cluster orientations and angular momenta, the number of 
particles per cluster should not drop below a minimum of 
a few hundred of particles. 

We search for bound systems using a friends-of-friends 
algorithm. In our ACDM model the virialization over- 
density p//5mcan at z = is ~ 330 ( Kitayama fc Sutc 
|1996| ) which corresponds (for spherical isothermal sys- 
tems) to a linking length of 0.17 times the mean inter 
particle distance. These bound systems are our clusters of 
galaxies, i.e. our cluster sized halos. The 3000 most mas- 
sive ones are the constituents of our mock cluster sam- 
ple. The mean distance between clusters in this sample is 
34.7 ft^Mpc, which is in agreement with the mean dis 



tance of clusters in the Reflex cluster survey (Bohringer 
et al. 1998). The most massive cluster has a mass of 



2.3 x 10 it> h- 1 M Q , resolved with 3700 particles, the light 



h ^ Ivxr. . . i r 

est still has 1.4 x lO 14 ft _1 M which corresponds to 224 
particles. 

The two-point correlation function £(r) of the cluster 
sized halos in redshift space shows the expected behav- 
ior (Fig. Q) compatible with the correlation function of 

~ (|2000D- 



the Reflex clusters determined in Collins et al 



The correlation length is 16 ft _1 Mpc, quantifying the am- 
plitude of the correlation function. At scales of about 
70 ft _1 Mpc £(r) becomes negative. 

The principal axes of a cluster are determined by 
means of the eigenvectors of the inertial tensor using all 
halo mass points. Assuming an ellipsoidal shape these vec- 
tors can be uniquely transformed into the three axes of the 
underlying mass ellipsoid. Subsequently the three axes are 
assumed to be ordered by size as a > b > c . In addition 
to the mass ellipsoid, the angular momentum vector L and 
the spin parameter A can be utilized to characterize the 
dynamical state of a cluster sized halo. L is calculated by 
summing up the angular momenta of all mass particles 
of the respective halo. The spin parameter is defined by 
A = uj/uj SU p, where to denotes the actual angular velocity 



and i 



'sup 



the angular velocity needed for the system to be 



rotationally supported against gravity. This ratio can be 
expressed as 



A = L/(G 1 / 2 M 3 / 2 i? 1 / 2 ), 



(1) 



where L is the absolute value of the angular momentum L, 
M the tot al mass of the system and R its radius, respec- 
tively (e.g. PadmanabhaE 1993 ). The radius R is taken as 
the radius of a sphere with the same volume as the halo 
(estimated on a fine grid). 
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Fig. 1. The two-point correlation function of cluster sized 
halos in redshift space (solid line). The dashed line has 
a slope of —1.8. The data points are taken from Collins 

(pooop . 



et al 



3. Mark correlation functions 

Mark correlation functions (MCFs) have been used for 
the characterization of marked points sets for quite some 



time (see 3toyan & Stoyan 1994 for an instructive intro- 
duction and overview). Bcisbart & Kcrschcr (2000) intro- 
duced MCFs to cosmology and used them to quantify 
the luminosity and morphology dependent clustering in 



the o b served galaxy d i stribu tion (see also Szapudi et al 



200C ) . Gottlober et al. ( 2002 ) discussed the correlations of 
galaxy sized ha los depending on the m erging history with 
MCFs (see also Gottlober et al. ( 2001 ) ) . In a recent review 
Bcisbart et al. ( 2002 ) show other physical applications of 
the MCFs and also how the MCFs can be calculated from 
stochastic models. To make this article self contained, we 
repeat briefly the basics of MCFs and then extend the 
formalism to allow for vector marks: 

Consider the set of N points {x^}^ and attach a mark 
nii to each point x,; G H. resulting in the marked point set 
{(xj, mi)}f =l flStoyan|[l984l ; |Stoyan fc Stoyan||l994| ). In the 



following we use a variety of marks: the mass, the absolute 
value of the angular momentum, the spin parameter A, and 
as vector valued marks, the principal axes of the second 
moment of the mass distribution and the angular momen- 
tum. Let p be the mean number density of the points in 
space and p M (m) the probability density of the mark dis- 
tribution. The mean mark is then m — j dm p M (m)m, 
and the variance is V[m] = J dm p M {m)(m — m) 2 . We 
assume that the joint probability p (x, m) of finding 
a point at position x with mark M, splits into a space- 
independent mark probability and the constant mean den- 



sity: p (m) x p. 

The spatial-mark product-density 

^((xi.miJ.Cxa.nfc)) dV 1 dm 1 dV 2 dm 2 , (2) 

is the joint probability of finding a point at xi with the 
mark mi and another point at x 2 with the mark m 2 . We 
obtain the spatial product density /9 2 (xi, x 2 ) and the two- 
point correlation function £(r) by marginalizing over the 
marks: 



P 2 (l + £(r))=/>a(xi,x 2 ) 



dmi / dm 2 pf M ((xx,mi), (x 2 ,m 2 )), (3) 



with £(r) only depending on the distance r = |xi — x 2 | of 
the points in an homogeneous and isotropic point set. 
We define the conditional mark density: 



X 2 (mi,m 2 |xi,x 2 ) 
(pi 



p 2 (xi,x 2 ) 







for /9 2 (x!,x 2 ) ^ 0, 
else . 



(4) 



For a stationary and isotropic point distribution, 
•M 2 (toi, m 2 |r)dmidm 2 is the probability of finding the 
marks mj and m 2 of two galaxies located at Xj and x 2 , un- 
der the condition that they are separated by r — |xx — x 2 |. 
Now the full mark product-density can be written as 



9-2 



>M ((xi,mi),(x 2 ,m 2 )) = 



= X 2 (mi,m 2 |xi,x 2 ) p 2 (xi,x 2 ). (5) 



If there is no mark-segregation A4 2 (mi,m 2 \r) is indepen- 
dent from r, and A4 2 (mi, m 2 |r) = p M (mi)p M (m 2 ). 

Starting from these definitions, especially using the 
conditional mark density M. 2 (rn\,m 2 \r) 1 one may con- 
struct several mark-correlation functions sensitive to dif- 



ferent aspects of mark-segregation ( Bcisbart fc Kerscher 
[2000D . The basic idea was to consider weighted correla- 
tion functions conditional that two points can be found 
at a distance r. For a non-negative weighting function 
f(mi , m 2 ) we define the average over pairs with separation 
r: 



(f) P (r) 



J dmi J dm 2 f(mi 



m 2 ) M 2 (mi,m 2 \r). (6) 



(/)p ( r ) i s the expectation value of the weighting function 
/ (depending only on the marks), under the condition that 
we find a galaxy-pair with separation r. For a suitably de- 
fined integration measure, Eq. (^) is also applicable to 
discrete marks. The definition (|J) is very flexible, and al- 
lows us to investigate the correlations both of scalar and 
vector valued marks. 

To calculate such mark correlation functions from the 
cluster data we use an estimator taking into account the 
periodic boundaries of the simulation box. We obtain vir- 
tually identi cal results for the estimator without boundary 
corrections ( Beisbart fc Kerschei 2000 ). 
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3.1. Correlations of scalar marks 

For scalar marks the follo wing mark correlation func tions 
have proven to be useful ( Beisbart fc Kerschcr 2000 ): 

— The simplest weight to be used is the mean mark, 
_ (mi + m 2 ) p (r) 



2 m 



(7) 



that quantifies the deviation of the mean mark of 
pairs with separation r from the overall mean mark 
m. A k m > 1 indicates mark segregation for point 
pairs with a separation r, specifically their mean mark 
is larger than the overall mark average. 

Accordingly, higher moments of the marks may be 
used to quantify mark segregation. The mark variance 
var(r) 



var{r) = ({m x - (m x ) p (r)) ) (r) 



(8) 

gives information about the fluctuations of the marks 
of points which are separated by a distance r to other 
members of the set. A var(r)/V[m] larger than one 
indicates a substantially increased scattering of the 
marks compared to th e overall vari ance. 
— The mark covariance ( Cressie 1991 ) is 

cov(r) ee ((mi - (mi)p (r))(m 2 - (m 2 ) P (r))) p (r) 

= (mim 2 ) p (r) - (mi) p (r) (m 2 ) P (r). (9) 

Mark segregation can be detected by looking whether 
cov(r) differs from zero. A cov(r) larger than zero, e.g., 
indicates that points with separation r tend to have 
similar marks. 

3.2. Correlations of vector marks 

To study the alignment effects between the clusters we at- 
tach to each cluster the orientation as a vector mark 1 with 
|1| = 1. The orientation is cither given by the direction of 
the major half-axis of the mass ellipsoid, or the direction 
of the angular momentum. The distance vector between 
two clusters is r, and the normalized direction is f = r/r. 
We consider the following MCFs: 



^(r) = (|l 1 .l 2 |) P (r), 

^(r) = i(|l 1 .f| + |l 2 -f|) P (r), 



(10) 



where "•" denotes the scalar product. Due to the spatial 
symmetry of the mass ellipsoids we use the absolute values 
of these scalar products. 

— A(r) quantifies the direct Alignment of the vectors li 
and 1 2 (see also [Stoyan k, Stoyaii 1994 and their k^). 



A(r) is proportional to the cosine of the angle between 
li and 1 2 . 

J-(r) quantifies the filamentary alignment of the vec- 
tors li and 1 2 with the line connecting both clusters. 
J-(r) is proportional to the cosine of the angle between 
li and the direction vector f connecting the points. 



In three dimensions no mark segregation implies A{r) = 
T(r) = 0.5, in two dimensions A(r) — T(r) — 2/tt. 
[Beisbart et al. (2002) provide some further explanations 
and discuss in which sense A(r), J-(r) provide a complete 
characterization of the vector correlations. 

4. Correlations in the orientation of clusters' 
shape 

As described in Sect. || we determine the mass ellipsoid 
for each of the cluster sized halos. The distribution of the 
intrinsic shapes of clusters in our sample is illustrated in 
Fig. ||. There we plot the ratios between the lengths of 
the three principal axes a > b > c: the upper panel shows 
the ratios a/b versus a/c, whereas the lower panel shows 
b/c versus a/c. For prolate rotational ellipsoids (a > b, 
b = c) all points would coincide with the dashed diagonal 
on the upper panel. In the case of oblate rotational ellip- 
soids (a = b, a > c) all data points would coincide with 
the dashed diagonal plotted in the lower panel. Clearly 
the majority of the data points tend to be concentrated 
rather along the upper panels' diagonal than along the 
lower panels' diagonal. Therefore, the clusters in our sim- 
ulations typically have prolate shape, in good agreement 
with observations (e.g. Cooray 2000) and other CDM sim- 
ulation, cp. (Cole & Lacey 1996). 

In the following we will use the direction of the ma- 
jor axis of the mass ellipsoid (the "a-axis") as a vector 
mark 1. Fig. ^| shows the mark correlation functions A(r) 
and F(r). The increased A(r) towards smaller scales in- 
dicates that pairs of clusters with a distance smaller than 
30 /i _1 Mpc prefer a parallel orientation of their orienta- 
tion axes li and 1 2 . This deviation of ~ 10% from the 
signal of a purely random orientation is clearly outside 
the random fluctuations shown by the shaded region. We 
determine these random fluctuations by repeatedly and 
randomly redistributing the orientation axes among the 
clusters. The positions of the clusters remain fixed. The 
filamentary alignment of the orientation of the clusters li 
and 1 2 towards the connecting vector f as quantified by 
F(r) is significantly stronger (up to 20% deviation from 
the signal of purely random orientation on small scales). 
Remarkably this signal extends out to 100 /i _1 Mpc. In a 
qualitative picture this may be explained by a large num- 
ber of clusters elongated in the direction of the filaments 
of the large scale structures. Filaments are prominent fea- 



tures found in the observed galaxy distribution (Huchra 



it al 



Shandarin 



1990) as well as in A-body simulations (Melott & 
199C), often with an extent up to 100 h^ 1 Mpc 



The inspection of density field of the simulation utilized 
here confirms this view (see Fig. ||). 

Up to now we used the orientation of the mass ellipsoid 
in three dimensional space. Most observations, however, 
only provide the projected galaxy number density or X- 
ray intensity. In the following we will investigate whether 
the clear signal found in Fig. |§| will be reduced by pro- 
jection. In close analogy to our three dimensional analysis 
we determine the orientation of the mass-ellipse of a halo 
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Fig. 2. Ratios of the length of the clusters principal axes 
(a > b > c). Upper panel: a/b versus a/c, for prolate ro- 
tational ellipsoids (a > b, b — c) all points would coincide 
with the dashed diagonal; lower panel: b/c versus a/c, for 
oblate rotational ellipsoids (a — b, a > c) all data points 
would coincide with the dashed diagonal. 



0.540 r 
0.530 r 




20 40 60 80 100 

r [ Mpc / h ] 




20 40 60 80 100 

r [ Mpc / h ] 



Fig. 3. Mark correlations using the 3D orientation of the 
dark matter halo, specified by the direction of the major 
axis 1 of the mass ellipsoid, as vector mark. The shaded 
area is obtained by randomizing the orientation among 
the clusters. 



from the mass-density orthogonally projected onto a side 
of the simulation box. We repeat our analysis with the 
major axis of the mass ellipse in two-dimensions as a two- 
dimensional vector mark. Also the normalized direction f 
is given in the plane. However, for the radial distance r 
we use the three dimensional separation of the clusters. 
With this setting, we mimic the observational constraints, 
e.g. in the Reflex cluster survey. In this projected sam- 
ple the signal in A(r) is only visible on scales smaller than 
10 ft. _1 Mpc (Fig. ||). The amplitude of the deviation from 
random orientation (in two dimension 2/tt) is reduced, as 
well. But still, the filamentary alignment quantified with 
F{r) indicates a strong alignment of the orientation of the 
clusters with the connecting vector (Fig. [|). Although the 
amplitude is slightly reduced, this alignment is still visible 
out to scales of 100 /i _1 Mpc, as in the three dimensional 
analysis. 

The question arises, why do we find two different 
regimes of the vector correlations, using either A{r) or 



3~(r). The direct alignment, quantified with A{r) ex- 
tends out to 30 h~ 1 Mpc, whereas the filamentary align- 
ment (.F(r)) is visible out to 100 h~ 1 Mpc. A glimpse 
on the real density distribution in Fig. [5] sheds light 
on this topic. Clearly one can see filamentary structures 
with sizes extending up to 100 h~ 1 Mpc, for example the 
two pronounced filaments which end up in the knot at 
x w 275 ft,~ 1 Mpc and y w 75 ft, _1 Mpc. Generically, such 
filaments are not straight but crinkled and show a lot of 
small branches. This implies that the coherence of the an- 
gles between the major axes of clusters harbored in these 
filamentary structures is lost soon. Consequently the pos- 
itive correlation seen with A(r) is confined to compara- 
tively small scales < 30 /i _1 Mpcs (in the case of projected 
data even to < 10 /i _1 Mpc). On the other hand, the sig- 
nal exploring the angle between orientation and connect- 
ing vector is not affected by the small scale disorder. Thus 
J-(r) is tracing the filamentary alignment out to scales of 
100 ft- x Mpc. 
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Fig. 5. Projected density field of a ~ 7.5 /i _1 Mpc slice 
through the simulation box. The density is calculated on 
a 300 3 grid which corresponds to a smoothing length of 
~ 1.7 /i Mpc. The maximum value of projected density 
is w 1.5 x 10 4 i o//9mcan. Gray scales are chosen according 
to log of density. 



Fig. 4. Mark correlations with the normalized major axis 
|1| = 1 of the projected cluster-density (2D orientation) as 
vector mark. The shaded area is obtained by randomizing 
the orientation among the clusters. 



5. Angular momentum, mass and spin correlations 

In a simulation we have access to the position and the ve- 
locity of the particles. Hence, it is easy to determine the 
angular momenta of our cluster sized halos. In this section 
we use the direction of the angular momenta of galaxy 
clusters as vector marks and compare their correlation 
properties with the alignment seen in the orientation of 
the mass distribution (Sect. [l]). With scalar MCFs, using 
the absolute value of the angular momentum, the clus- 
ter mass and the spin parameter as scalar mark, we will 
supplement the foregoing result. 

It is well known that there exists a correlation between 
the axes of the mass ellipsoid and the direction of the 
angular momentum in gravitationally bound N-body sys- 
tems ( Binncy fc Trcmainc| 1987). Fig. ^| shows that the 
major axis of the mass ellipsoid tends to be perpendicular 
to the angular momentum. Compared to the expectation 
from a purely random distribution, the minor axis favors 
smaller angles with the angular momentum. Hence, using 
the direction of the angular momentum as vector mark, 
we are expecting to see effects comparable to the shape 



_ 0.8 




Fig. 6. Cumulative plot of the number of clusters and the 
cosine between the direction of the angular momentum 
and the major, medium and minor cluster axis, respec- 
tively. 

orientations. The vector correlations of the directions of 
the angular momentum are shown in Fig. [7| There is no 
clear signal for a direct alignment of the angular momenta 
as traced by A(r). However, the correlation F(r) < 0.5 is 
smaller than for random directions (shaded area), this in- 
dicates that the angular momenta are preferably oriented 
perpendicular to the connecting vector f on scales up to 
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Fig. 7. Mark correlation functions using the normalized 
angular momentum 1 with |1| = 1 as vector mark. The 
shaded area is obtained by randomizing the orientation 
among the clusters. 
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Fig. 8. Mark correlation functions with the absolute value 
of the angular momentum of the cluster as scalar mark. 
The shaded area is obtained by randomizing the mark 
among the clusters. 



30 /i _1 Mpc. The amplitude of the deviation from the ran- 
dom orientation reaches approximately 5% and is clearly 
outside the fluctuations. With Fig. ^ in mind, the latter 
result was expected from the correlations of the orienta- 
tions shown in Fig. For a pair of clusters residing in 
a filamentary structure the connecting vector is oriented 
along this filament. The major axes of the orientations 
tend to align with the filaments, and therefore the angular 
momenta tend to stand perpendicular to the filaments. A 
random orientation of the angular momenta in the planes 
perpendicular to the filaments provides a simple explana- 
tion for the absence of any signal in A(r) , still compatible 
with the strong signal in J-{r). 

Besides the correlations in the directions of the angu- 
lar momentum vectors, we are also interested in the cor- 
relations of their magnitudes. We investigate these corre- 
lations with the MCFs discussed in Sect. |3.l[ using the 
magnitude of the angular momentum as a scalar mark. 
In Fig. H the increased k m (r) shows that pairs of clusters 
with separations < 50 /i _1 Mpc tend to have higher mean 
angular momentum compared to the overall mean angular 



momentum. The positive covariance up to < 15 h~ Mpc 
indicates that both members of close pairs tend to have 
similar angular momentum. Hence, inspecting both fc m 
and the covariance, we see that close clusters tend to carry 
a similar amount of angular momentum, larger than the 
overall mean angular momentum. 

There are two possibilities to explain this behavior. On 
the one hand the mean mass of a close cluster pair could 
be enhanced, and so the absolute value of the angular 
momentum would grow simply due to the fact that bigger 



clusters are under consideration. Indeed. Gottlobcr et al 



(2002) found such an increase in the mean mass in close 
pairs (in their case traced by the maximum circular veloc- 
ity of the halo). Consistently Beisbart & Kerscher (2000) 



report an enhanced luminosity for close galaxy pairs. On 
the other hand, the rotational support of close cluster pairs 
could be higher, meaning that the spin parameter A is en- 
hanced, cp. Eq. ([!]). To illustrate this further, we investi- 
gate the MCFs with the mass and the spin parameter as 
marks: 

Mass: In Fig. ^ we show the MCFs of the cluster distri- 
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Fig. 9. Mark correlation functions with the mass of the Fig. 10. Mark correlation function with the spin parame- 
cluster halo as scalar mark. The shaded area is obtained ter A of the cluster halo as scalar mark. The shaded area 
by randomizing the mark among the cluster halos. is obtained by randomizing the mark among the clusters. 



bution with the total mass as scalar mark. The increased 
k m (r) indicates that close pairs of clusters tend to have 
higher mean masses than the overall mean mass m on 
scales out to 50 h Mpc. The signal shows a deviation 
up to 10% from a purely random distribution. The condi- 
tional covariance cov(r) shows only a weak positive signal, 
confined to scales below 15 /i _1 Mpc, indicating that only 
close pairs tend to have similar masses. 
Spin: In Fig. |o| we show the MCFs of the cluster dis- 
tribution with the spin parameter A as scalar mark. The 
increased k m (r) indicates that neighboring pairs of clus- 
ters tend to have higher spin parameter A compared to 
the overall mean spin A. The signal seen in Fig. |l^ shows 
an enlarged spin parameter of clusters on scales out to 
50 ft, Mpc, the same range as for the mass. The deviation 
from random distribution is ~ 15%. Hence, not only an 
increased mass accretion but also tidal interactions out to 
fairly large scales seem to characterize the regions around 
cluster halos. The conditional covariance cov(r) shows no 
signal. This can be explained by looking at the mark vari- 
ance var(r). The signal at small distances deviates about 
50% from random, indicating large scattering of the spin 



parameter. This large dispersion also might be the rea- 
son why previous studies did not find an environmental 
dependence of the amount of the spin parameter. A suffi- 
ciently large simulation and also a large number of clusters 
are needed to get a clear signal in the MCFs for the spin 
parameter. 

Putting together all the results shown in Figs. |s|-[io|, 
it turns out that the increased mean angular momentum 
is not only caused by an enhanced mass, but also by an 
enlarged spin parameter of neighboring pairs of clusters. 
From the present analysis we can not draw any firm con- 
clusions, whether this phenomenon is dominated by tidal 
interaction or merging processes. There are compelling 
observational hints (see e.g. Plionis & Basilakoe 2002; 
[Sclmccker et al. 2001) that clusters in high-density envi- 
ronments show indications of dynamical disruption. Thus 
it seems to be likely that the enhancement of the angular 
momentum of close neighbors is caused by a substantial 
mass accretion (i.e. merging) . The concordant behavior of 
the mass and the spin MCF support this interpretation, 
at least on scales below ~ 15 h~ 1 Mpc. However, tidal in- 
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teractions could be the cause for the correlations seen on 
large scales. 

6. Summary 

Whether there exist correlations in the orientations of 
galaxies or galaxy clusters has been discussed for a long 
time. Binggeli ( 1982] ) reported a significant alignment of 



the observed galaxy clusters out to 50 h 1 Mpc. Struble & 



Peebles ( 1985 , 1986 ) clai med that thi s effect is small and 
prone to systematics and Ulmcr ct al. ( 1989 ) find no indi- 
cation in their investigation. Subsequently, several authors 
found, sometimes only weak, signs of alignments in the 



galaxy and galaxy cluster distribution (see e .g. Djorgovski 



1986: Lambas et al. 1988: Fuller et al. 199E; Heavens et al 



200C). As a novel statistical method we have used the 



mark correlation functions (MCFs) to quantify the align- 
ment of cluster sized halos, extracted from a large scale 
simulation based on a ACDM cosmology. Our sample with 
3000 cluster sized halos is bigger than the currently avail- 
able samples of galaxy clusters. The unambiguous signal 
we obtain benefits from the large statistics in our simula- 
tion. 

Using two different weighting functions in the con- 
struction of the MCFs we investigate the direct align- 
ment and the filamentary alignment. First we use the 
major axis of the mass ellipsoid as our direction marker. 
The clear signal from the direct alignment A{r) extends 
out ~ 30 h~ 1 Mpc. For the filamentary alignment T(r) 
we find deviations from isotropy up to ~ 100 /i _1 Mpc. 
Considering the projected mass distribution, the signal 
from the direct alignment A(r) already vanishes at a scale 
of ~ 10 h Mpc. However, we find a filamentary align- 
ment T(r) out to scales of ~ 100 /i _1 Mpc, even for the 
projected data. This scale is very similar to the size of 
the large scale filaments seen in our simulation. We think 
that the function T(r) is a powerful tool for exploring 
large scale alignment effects also in observational data. 



Franx et al. (1991) showed that the angular momentum 
of an ellipsoidal system tends to align with the minor axis 
of this system. We confirm this behavior in our simulation. 
With the angular momentum as vector mark, !F(r) shows 
the expected filamentary correlations: the angular momen- 
tum tends to be perpendicular to the connecting line, i.e. 
the filament, up to separations of ~ 40 h~ 1 Mpc. However, 
we obtain no signal for the direct alignment A(r). This is 
in concordance with the perception that the angular mo- 
menta are randomly oriented in the planes perpendicular 
to the filaments. 

With the scalar MCFs k m (r) and cov(r) we have inves- 
tigated the correlations in the absolute value of the angu- 
lar momentum. Close pairs of clusters tend to have similar 
and also higher absolute values of the angular momentum 
compared to the global average. A clear signal can be de- 
tected up to ~ 50 /i _1 Mpc. A further analysis of the mass 
and spin parameter distribution of the clusters with the 
MCFs has shown that this enhancement of the absolute 
value of the angular momentum is caused by an enhanced 



mass of close pairs of clusters as well as by the stronger ro- 
tational support of them. This behavior should be caused 
by the combined action of large-scale tidal fields and the 
hierarchical merging of progenitor structures and mass in- 
flow onto the cluster. Since this mass growth follows the 
large scale filaments, tidal interactions and merger events 
are tightly connected. The mark correlation function with 
scalar and vector marks deliver quantitative measures of 
these effects. 
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